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Abstract 

Plane-wave representations are used to formulate the exact solutions to frequency-domain and 
time-domain sources illuminating a magnetodielectric slab with complex permittivity e(u>) and 
permeability In the special case of a line source at z = 0 a distance d < L in front of an 

L wide lossless double negative (DNG) slab with k(uj 0 ) = e(w 0 )/e 0 = //j, 0 = —1, the single¬ 

frequency (cj 0 ) solution exhibits not only “perfectly focused” fields for z > 2L but also divergent 
infinite fields in the region 2d < z < 2L. In contrast, the solution to the same lossless k(wq) = — 1 
DNG slab illuminated by a sinusoidal wave that begins at some initial time t = 0 (and thus has 
a nonzero bandwidth, unlike the single-frequency excitation that begins at t = —oo) is proven to 
have imperfectly focused fields and convergent finite fields everywhere for all finite time t. The 
proof hinges on the variation of k(uj) about co = cu 0 having a lower bound imposed by causality and 
energy conservation. The minimum time found to produce a given resolution is proportional to 
the estimate obtained by [Gomez-Santos, Phys. Rev. Lett., 90 , 077401 (2003)]. Only as t —> oo do 
the fields become perfectly focused in the region z > 2L and divergent in the region 2d < z < 2L. 
These theoretical results, which are confirmed by numerical examples, imply that divergent fields 
of the single-frequency solution are not caused by an inherent inconsistency in assuming an ideal 
lossless k(u 0 ) = — 1 DNG material, but are the result of the continuous single-frequency wave 
(which contains infinite energy) building up infinite reactive fields during the infinite duration of 
time from t = —oo to the present time t that the single-frequency excitation has been applied. An 
analogous situation occurs at the resonant frequencies of a lossless cavity. A single-frequency (zero 
bandwidth) source inside the cavity produces infinite fields at a resonant frequency, whereas the 
same source turned on at time t = 0 (so that it has a nonzero bandwidth) produces finite fields. 

PACS numbers: 41.20.Jb, 42.25.Bs, 42.25.Fx, 42.30.Kq 
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I. INTRODUCTION 


The main purpose of this paper is to explain and resolve a number of the peculiarities 
and apparent paradoxes (such as perfectly reproduced source fields as well as divergent fields 
to the right of the slab) exhibited by the solution to a single-frequency (zero bandwidth) 
sinusoidal source illuminating a lossless magnetodielectric slab with relative permittivity and 
permeability equal to negative one [1]—[6]. This is accomplished by determining the solution 
for a time-domain (nonzero bandwidth) source produced by turning on the same sinusoidal 
source at some finite initial time [4], [6], [7], for example, at 7 = 0, and assuming the lossless 
slab has the slowest possible frequency variation in relative permittivity and permeability 
(about the value of negative one) allowed by causality and energy conservation. To clearly 
reveal the peculiarities and apparent paradoxes in the single-frequency solution, however, we 
begin by deriving a rigorous plane-wave solution to the single-frequency source illuminating a 
general lossless or lossy magneto dielectric slab with arbitrary permittivity and permeability. 


II. FREQUENCY-DOMAIN SOLUTION 


The boundary value problem of a time-harmonic tv > 0) source illuminating an 

infinite magnetodielectric slab can be solved simply and rigorously in terms of plane-wave 
representations [8]—[12]. For example, the plane-wave solution for the x component of the 
electric field E x (x,z) of a transverse electric (TE) (E y = 0, H y ^ 0) line source with no 
variation in the y direction located a distance d in front of a slab (infinite in the x and y 
directions and normal to z ) with width L, complex permittivity e = e' + ie", and complex 
permeability fi = y! + ifi" is given by (see Figure 1) 


where 


E x (x,z) 


+oo 



7o = ( k o 
7 = (k 2 


T 0 (h)e ^° 2 + R 0 (h)e-^ z 

, 0 < z < d 


T s (h)e ijz + R s {h)e~^ z 

, d < z < d + L 

(i) 

T(h)e^° z 

, d + L < z 


-h 2 )K kl = u 2 fi 0 e 0 


(2a) 

— h 2 )^ , k 2 =u 2 ye 


(2b) 


and e 0 and /i 0 are the permittivity and permeability of the free space in which the slab is 
assumed located. For passive materials e" > 0 and ji" > 0. The square root in the definition 
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( 2 a) of 7 0 is chosen positive real or positive imaginary depending upon whether h 2 < k% or 
h 2 > hg, respectively. The sign of the square root in the definition in (2b) of 7 is chosen 
to keep the imaginary part of 7 positive. If k 2 is real and h 2 < k 2 , then 7 is real and the 
sign of 7 is found by inserting a small loss, choosing the imaginary part of 7 positive, and 
letting the loss approach zero. This procedure leads to a positive real 7 if e and 11 are both 
positive real and a negative real 7 if e and y are both negative real. (The signs of the square 
roots can also be determined from the requirement that the energy flow in the incident and 
transmitted propagating plane waves and the held decay in the incident and transmitted 
evanescent plane waves be away from the source.) 



The plane-wave spectrum T 0 (h ) of the fields to the right (z > 0) of the line source (incident 
fields) is assumed given. For example, assume the TE line source is a two-dimensional y 
directed magnetic line current (magnetization). Then T 0 (h) is independent of h and can be 
written as 

Mh) = (3) 

K 0 

where E 0 is a constant with electric held units and the constant A ; 0 is inserted into the 
denominator of (3) to ensure the dimension of E x (x,z) is explicitly that of an electric held. 
If at some frequency cn 0 , the constitutive parameters fj,(u 0 )/fj, Q = e(u 0 )/e 0 = —1, it will be 
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shown that the slab-induced helds of the y directed magnetic-current line source diverge to 
infinite values in certain regions. Moreover, this infinite divergence is not peculiar to that 
particular source. 

The reflected spectrum Rq{K) to the left of the slab, the transmitted and reflected spectra, 
T s {h) and R s (h ), within the slab, and the transmitted spectrum T(h) to the right of the 
slab are obtained by equating the tangential components of the electric and magnetic helds 
across the interfaces of the slab at z — d and z = d + L. Specifically, 


T a {h) = 

R s (h) = 
Ro(h) = e 


T(h) = T 0 (h)T TE (h) 

Tity ( i , £o7\ i( 7o _ 7 )(d+i) 

2 l W 

T ( h ) f 1 _ e h7o+7 )(d+L) 

2 V 7 tf) 

ho d ( Tse hd + Rge -ind _ Toe n 0 dj 


where the TE transmission coefficient is given as [2] 

4p-n 0 L 

T TE (h) = 


( 2 + f2o + 5ol) e -i 7 L + (2 - nsL 
V e o7 n 0 J v e o7 


€ (A A g*7^ 
e 7 0 / 


(4a) 

(4b) 

(4c) 

(4d) 
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A. Lossless —1 double-negative slab 


In the case of the “perfectly focusing” slab [1], [2], e/e 0 = n/n 0 = —1 at some frequency 
uJq and we have 7 = —q 0 if h 2 < k 2 — kjj 0 = caQ/i 0 e 0 and 7 = 7 0 if h 2 > k 2 ]0 . Then 


T„(h) = e-^ L 


( 6 ) 


and (4) become 



T(h) — 

T 0 (h)e~ 

-i7 0 2L 




T.(h) = | 

\ T 0 (/r)e i7 ° 2d 

, 7 = 

= —7o , 

h 2 

< 

k 2 

^00 

1° 

, 7 = 

= 7o , 

h 2 

> 

otc 

0 

R s (h) = < 

f° 

\ _ . 

, 7 = 

= -7o > 

h 2 

< 

k 2 

^00 

y T 0 (h)e^° 2d 

: , 7 - 

= 7o , 

h 2 

> 

O 
^ O 


Ro(h) — 0 


(7a) 

(7b) 

(7c) 

(7d) 
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so that E x (x,z) in (1) can be written as 


E x {x,z) 


^ J dhT 0 (h)e ihx l 


e ilo z 

e il 0 (2d-z) 

e *7oO-2£) 


, 0 < z < d 
, d < z < d + L 
, d + L < z . 


( 8 ) 


Since the incident helds are equal to the helds to the right of the source in free space, that 

is 

+oo 

E'r{x,z) = — [ T 0 (h)e^ hx+ ^dh, z> 0 (9) 

2n J 

—oo 

by referring to (3), the equations in (8) can be re-expressed as 


' 

E x (x,z ) = < 


E™(x, z) , 0 < z < d 

E™ c (x, 2d — z) , d, < z < 2d 

infinite divergent held , 2d < z < 2L 
E™(x, z-2L ) ,2 L< z 


( 10 ) 


where from herein out it is assumed that d < L. In other words, the held to the right of 
the source and to the left of the slab is just the incident held of the source in free space. 
The held to the right of the front face of the slab and to the left of z = 2d is the image 
of the source held. The held between z = 2d and z = 2L diverges to infinite values. Most 
importantly, and quite remarkably, the held in free space to the right of z = 2L is just 
the incident held translated to the right a distance equal to twice the width of the slab. 
The phase and magnitude changes of the propagating and evanescent plane waves in the 
free-space regions between 0 < z < d and d + L < z < 2L are canceled by opposite phase 
and magnitude changes in the —1 double negative (DNG) slab. This perfect replication of 
the free-space source helds for z > 0 in the free-space region z > 2L to the right of the slab 
is sometimes referred to as “perfect focusing” [1], [2], 


B. Lossy —1 double-negative slab 

The solution in (8)—(10) is so unusual that it warrants further investigation. It is a solution 
that assumes the loss in the slab (and the surrounding space) is exactly zero. It may be 
physically more appealing to insert a small loss into the slab [17], [18], [23] and determine the 
solution for the helds as the loss approaches zero. This solution can be expressed rigorously 
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from (1) in terms of two limits: the infinite limit for the evanescent spectrum and the limit 
as the loss in the slab approaches zero. For simplicity, let the relative loss at the frequency 
cu 0 in /x(cu 0 ) and e(cu 0 ) of the slab be equal and denoted by 


^ (^o) — h O^oV/h) ~ 6 ( UJ o)/ e o > 0- 


( 11 ) 


Then we can write from (1) 


Et 


+H 


— lim lim 

2 n 6"^0 oo 


dhe 


ihx 


-H 


T 0 (h)e iloZ + Ro(h)e~ iloZ , 


T s (h)e ilz + R s {h)e~^ z 
T(h)e iloZ 


0 < z < d 
d < z < d + L 
d T L < z . 


( 12 ) 


One could argue that, in practice, the wavenumber \h\ of the evanescent spectrum incident 
upon the slab should always be truncated to a finite limit H 0 > /c 0 o = ^ox/ho^o because the 
evanescent spectrum for \h\ greater than some H 0 > k 00 will be lost in the noise. Then the 
lim h^oo flu would be replaced by merely , the lim^/^o could be brought under the 
integral sign, and the solution in (8)-(10) would be approached for H 0 k 00 . 

Still, one could ask what the solution becomes if, in principle, the evanescent spectrum 
is not truncated and the lim^/_» 0 is not brought under the integral sign in (12). In that 
case, we have e 7 0 /(e 0 7 ) = —1 + id"{ 1 + ^oo /1 To 1 2 ) + 0[(d") 2 ] ~ —1 + id" + 0[(d") 2 } for the 
evanescent spectrum if terms in A)q 0 /| 7 0 | 2 are neglected compared to unity and we find that 
Z rE (/r) in (5) for 5" -C 1 can be approximated by 


T 

TE 


(h) 


ho\ L 


e r '° 


^ e h 0 \L + e -h 0 \ L ’ 


h 2 \ 

u ^ n 00 . 


(13) 


This expression reveals that Z rE (/r)e _ bwhich comprises the evanescent integrand in (12) 
in the region z > d + L, rapidly decreases toward zero as | 7 0 | grows larger than Ta, where 
Ta is given implicitly by 


6" « e“ r ^ /2 , z > d + L . (14a) 

Choosing the minimum value of z = d + L, this expression for Ta becomes 

5" w e -rs(d+L)/2 ( 14b ) 

so that 

r 5 ~ — -j~~r In d" . (14c) 
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Moreover, X rB (/r) is given approximately by (6) for the entire propagating spectrum ( h 2 < 
/cq 0 ) with 6" < 1, and for the evanescent spectrum up to that is, 0 < | 7 0 | < IV Thus, 
for the evanescent spectrum in the domain /c 0 o < h < , where 


H5 ~\j{dh lnS ") +k ° 2 ° (15) 

(4) becomes 

T(h) « T 0 (/r)e 2|7o|i (16a) 

T s (/i) « —iT 0 (h)e < ' 2|7ol ” r,5 - )L (16b) 

R a (h) « T 0 (h)e- 2|7 o |d (16c) 

i? 0 (h) « -iT 0 (h)e 2 l 7 oK L - d )- r5 - L (16d) 


and the evanescent part of the held of (12) becomes 


E7’ S "-°(x,z) 


lim — / dhT 0 (h)e 
5"->o 2vr ,/ V 7 

&00<|/l|<#c5 




e -|7oh _ j e h 0 \(z+2L-2d)-r s L 
g |7 o (z—2d) — i ^^\'Iq\[2L—z)—T^L 
e l7 0 l( 2i —d 


, 0 < z < d 
, d < z < d + L 
, d + L < z 


(17) 

where Hs is given in (15). For z<A 0 (where A 0 = 27T//c 00 denotes the free-spaee wavelength), 
that is, z extremely close to the source, Hs may have to be increased to include all the 
significant evanescent waves. However, if 2 > 2 L, the value of z in (14a) can be chosen equal 
to its minimum value of 2 L in that region to obtain 


5" « e~ VsL , z>2L 


(18a) 


and Ts in (14c) can be replaced by 

T s « — In 5", z > 2L 

Ju 


so that Hs becomes (see also [17], [18], [23]) 



z > 2L. 


(18b) 


(19) 
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The propagating spectrum is the same as in (8) as 5" —> 0. Therefore, the propagating 
spectrum in (8) combines with the evanescent spectrum in (17) to yield 

bounded field , 0 < z < d — (L — d) 

infinite divergent field , d — (L — d) < z < d + (L — d) 

E x ~^°(x, z) = < bounded field , d + (L — d) < z < 2d (20) 

infinite divergent field , 2d < z < 2L 

, E™(x,z-2L) ,2 L<z 

where it is assumed in (20) that L/2 < d < L. The fields in (20) conform to those obtained 
from the solution in [13], [14]. 

Comparing (20) and (10) reveals that the fields in the region 0 < z < 2d differ depending 
upon whether the loss 5" in the slab is made to approach zero before [to get(10)] or after 
[to get(20)] the limit of the integration of the evanescent spectrum is allowed to approach 
infinity. However, in the important free-space region to the right of the slab (z > d + L), the 
nature of the fields is independent of the order in which the limit of the loss (approaching 
zero) and the limit of the integration of the evanescent spectrum (approaching infinity) is 
taken. In either case, “perfect focusing” of the source fields for z > 0 is attained in the limit 
as the loss approaches zero in the free-space region z > 2L to the right of the slab. 

Also, in either case, as the loss approaches zero, the field diverges to infinite values in the 
free-space region d + L < z < 2L that lies to the right of the slab. The fields throughout 
the free space to the right of the slab (z > d + L) may, at first sight, appear to violate 
the analyticity theorem [15, ch. V, sec. 4], [16, Theorem 2.2], which states (to quote [16, 
Theorem 2.2]), “If u [our E x field] is a two times continuously differentiable solution to 
the [homogeneous] Helmholtz equation in a domain T> [our region z > d + L\, then u is 
analytic.” [24] This theorem seems to imply that the fields in the free-space region to the 
right of the slab (z > d + L) should be analytic, whereas in fact they are analytic in the 
region z > 2L, but diverge to infinite values in the region d + L < z < 2L [5]. This apparent 
paradox is resolved if it is noted that the analyticity theorem requires that the function 
be a twice continuously differentiable solution to the homogeneous Helmholtz equation. It 
may be possible to weaken this condition to something less restrictive, but certainly the 
theorem does not apply to fields that diverge to infinite values in part of the region, namely 
(d + L < z < 2L) — as indeed our plane-wave solution demonstrates. [25] In particular, 
an electric field component in a free-space region need not be an analytic function of the 
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spatial coordinates throughout that region if the held is allowed to diverge to infinite values 
in part of that region. 

For a small but nonzero value of the loss parameter 8", we note from (IT) [before 8" — > 0 
in (17)] that the fields everywhere to the right of the source (z > 0) are bounded. Moreover, 
with a small but nonzero loss 8", the fields throughout the free-space region to the right of the 
slab ( z > d + L), and throughout the region between the source and the slab (0 < z < d), 
are analytic functions of complex x and z in complex neighborhoods of the real x and z 
coordinates. For a nonzero loss 8", equation (17) shows that the held in the free-space 
region z > d + L to the right of the slab is given approximately by 

+h s 

E 5 Jx,z)tt^- [ T 0 (h)e i[hx+ ^ z ~ 2L)] dh, z>d + L (21) 

2tt J 
-H S 


with Hg given in (15) for d + L < z < 2L and in (19) for z > 2 L. The transverse resolution 
just to the right of z = 2L can be found by integrating (21) for z = 2L with T 0 (h) for a 
magnetic-current line source inserted from (3) to get 


Ei(x,2L) 


E 0 sin (Hgx) 


IT k 


oo 


x 


( 22 ) 


which shows that the held of the line source is approximated by a sine function with a waist 
size proportional to 1 /H§. For two identical line sources separated along the x axis an equal 
distance D/2 from the origin, the x component of the electric held at z = 2L is 


sm[H s (x- D/2)] | sm[H s (x + D/2)]] /ooA 

x-D/2 + x + D/2 J ' 1 j 

Numerical computations of the held in (23) of the two sine functions show that the resolution 
Ax, defined as the separation distance such that the two peaks produced by the two sine 
functions are 3 dB in intensity above the central minimum, is given by [26] 


Ex’ D ( x , 2L) 


En 


nk, 


oo 


Ax 


1.53?r 1.53vr 



(24) 


The “resolution enhancement” R e , defined as the ratio of the resolution with H$ to the 
resolution with the propagating waves alone (Hg = k 00 ), is thus given simply as 


Re 


El 

koo 



(25) 
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so that 


(26) 


. 1.537T .76 

Ax = -—— = — A 0 . 


k 00 Re Rp 


For — Aq \a.8"/(2 , kL) S> 1, (25) reduces to 


R e 


A 


2nL 


° Ini" 


(27) 


an expression for resolution enhancement derived by Smith et al. [18]. The 5" in (25) needed 
to obtain a given resolution enhancement R e is 


S" « 


(28) 


For example, to obtain a resolution enhancement of R e = 5 in a one wavelength slab (L = 
Ao), the loss parameter 5" should have a value no larger than about 


8 " « e - 2?rv/ 24 


4.3 x 10 


-14 


(29) 


which is an extremely small loss. For R e = 2.5 resolution enhancement, the loss should be 
no larger than about 

. (30) 


8" w e 27 r V (2 - 5 ) 2 1 ps 5.6 x 10 7 


quite a small yet more realistic value [19]. The resolution formulas (24)-(26) are confirmed 
by the numerical examples in Section IV. 

As 8" —> 0 the increasingly large fields in the region d + L < z < 2L can be evaluated 
asymptotically from (21), (19), and (3) for the y directed magnetic-current line source to 
get (for z not too close to 2 L) 

E 0 1 


E S x (x,z) « 


-- . „ , T cos I — In 5" + tan 1 

nk 00 ^x 2 + (2L - z) 2 {8") 2 -R L \L 2 L-z 


x 


X 


(31) 


The asymptotic {5") 2 Z R decay in (31) agrees with that found for the fields in [13]. 


III. TIME-DOMAIN SOLUTION 

In principle, the frequency-domain plane-wave solution to the lossless infinitely long mag¬ 
netodielectric slab with e(u; 0 )/e 0 = y(u 0 )/y 0 = —1 has shown that such a slab reproduces 
the incident fields of a single-frequency (zero bandwidth) source within the free-space region 
z > 2L to the right of the slab. We have shown that this result holds regardless of whether 
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the loss in the slab material is set equal to zero throughout the formulation of the solution 
or the loss is chosen nonzero during the formulation of the solution and then allowed to ap¬ 
proach zero. In either case, however, this single-frequency plane-wave solution also predicts 
that the fields diverge to infinite values in certain regions (defined in ( 10 ) or ( 20 )) between 
z — d — (L — d) and z = 2L with d < L. 

In practice, there will always be some loss in the slab material, some inhomogeneities 
within the slab material, and some noise level that limits to a finite value the effective 
wavenumber of the source evanescent waves that reach the front face of the slab. Also, a 
realistic slab will always be limited to a finite length. These practical realities will reduce the 
fields everywhere to finite values. Nonetheless, it is somewhat disconcerting that the exact 
classical solution to a source in front of an ideal lossless infinitely long magnetodielectric 
slab with e(u 0 )/e 0 = /i(a; 0 )//x 0 = —1 has infinite held values in certain regions defined in 
(10) or (20). These infinities may lead one to question the possible existence even within 
classical physics of an ideal lossless material with e(u; 0 )/e 0 = /i(u; 0 )//i 0 = —1 [3], [5]. 

In this section we will illuminate the slab with a time-domain (nonzero bandwidth) si¬ 
nusoidal wave that turns on at a given finite time in the past (as would any signal in the 
laboratory) [4], [ 6 ], [7], unlike a single-frequency (zero bandwidth) sinusoid that begins in 
the remote past (t = — oo) and thus illuminates the slab for an infinite amount of time. The 
former time-domain sinusoid carries a finite amount of energy between the time it turned 
on and the present time t, whereas the latter single-frequency sinusoid carries an infinite 
amount of energy between the time it turned on (t = —oo) and the present time t. The solu¬ 
tion that we obtain in this section to the time-domain sinusoid that is turned on in the finite 
past reveals that at a finite present time t all fields in the lossless e(u 0 )/e 0 = fj,(u 0 )/fj, 0 = —1 
slab are finite in value, and only as the present time t —> oo do values of the fields approach 
infinity in certain regions (the same regions where the single-frequency fields diverge to in¬ 
finite values). In other words, it is not the lossless slab with e(u 0 )/e 0 = /x(cu 0 )//i 0 = —1 
that inherently leads to the infinitely large single-frequency fields, but the single-frequency 
continuous wave that illuminates the slab from t — — oo to the present time t and imparts 
an infinite amount of reactive energy in certain regions within and near the slab. The time- 
domain solution unequivocally explains the origin of the infinite fields encountered in the 
lossless single-frequency solution. 

An analogous situation occurs for a sinusoidal source inside a perfectly conducting (that 
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is, lossless) cavity. If a single-frequency source is placed inside the cavity, the fields are well- 
behaved except at the resonant frequencies of the cavity where the fields diverge to infinite 
values. One does not conclude from these infinite divergences that a lossless cavity cannot, 
in principle, exist. One simply acknowledges that the infinite energy in the continuous 
wave has led to infinite reactive fields in the lossless cavity at the resonant frequencies. If 
the single-frequency source is replaced by a time-domain sinusoidal source that begins at a 
finite time in the past, one finds that the fields inside the cavity remain finite for all finite 
present time t even at the resonant frequencies. Only as t —> oo do the values of the fields 
of the cavity at the resonant frequencies approach infinity. 

A. Time-domain solution to the right of the slab (z > d + L) 

The frequency-domain solution for the TE line source in the region z > d + L can be 
rewritten from (1) and (4a) as 

+OO 

E x (x, z) = — f T 0 (h)T TE (h)e^ hx+ ^ z) dh , z > d + L (32) 

2vr J 

— OO 

or, alternatively, with the change of integration variable u = h/uj 

+oo 

E x {x,z) = — / T(tt)T TE (uM)e iw( “ +c o z) d«, z>d + L (33) 

27 T J 

— OO 

where u( 0 = 7 0 (c ou) so that 

Co = (l/c 2 -u 2 )^, /i 0 e 0 = l/c 2 (34a) 

with c being the speed of light in free space and 

T (u) = uT 0 (uu). (34b) 

It is assumed that T(u) is independent of the frequency ui, for example, as it would be for 
the y directed magnetic-current line source given in (3), that is 

T (u) = uT 0 (uu) = ——^ = E 0 c. (35) 

fi'o 

Equation (33) is the single-frequency (ca > 0) solution in the region z > d + L to the 
right of the slab for a line source (at z — 0) with e~ lut time dependence that has existed 
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from t = — oo in the remote past. For a sinusoidal wave, cos(oo 0 t), that turned on at a finite 
time t = —to in the past and turns off at some future time t = +f 0 , the frequency-domain 
spectrum T w is given by 

+£o 

oo Q t) e lLVt dt = — 

0 ; 2 tt 

-to 

in which we assume to > 0 and oo 0 > 0. (Choosing the start and end times of the cosine 
wave equal to —to and +to, respectively, simplifies the time-domain analysis. At the end of 
the analysis, the start and end times will be changed to 0 and t, respectively.) 

The x component of the time-domain electric held E x (x,z,t) can now be found by mul¬ 
tiplying the integrand in (33) by the frequency spectrum Tb and taking the inverse Fourier 
transform with respect to oo. Using the fact that E x (x,z,t) is a real function allows one to 
integrate over only the positive frequencies in the inverse Fourier transform such that [12, 
sec. 5.3] 

E x (x, z, t ) = Re [E x (x, z, t)} (37) 



sin [(a; + oo 0 )t 0 ] sin [(a; - a> 0 )t 0 ] 


00 + C On 


00 — OOn 


(36) 


where 

+oo +oo 

E x (x, z,t) — J f T^T^T^uuy^+^-^duduo, z>d + L. (38) 

0 —oo 

is the analytic-signal time-domain electric held [12, sec. 5.3]. Note that as t 0 —> oo, the 
frequency spectrum T u —> [<5(a; + a; 0 ) +<5(ct; — a; 0 )]/2, where S(x) denotes the delta function, 
since 


1 sin[(uo ± w 0 )i 0 ] 1 r/ , \ 

— lim --- = -o{u±Ur ] ) 

2n t 0 ^oo oo EoOq 2 

Thus, for large t 0 (38) can be expressed approximately as 


(39) 


+oo +oo 


E x (x, z, t) 


2n 2 



sin [(a; — a; 0 )to] 


00 — 00n 


T(u)T te (uovje^+^-Qdu doo , z>d + L. (40) 


0 —oo 

For large to, the dominant contribution in (40) will come from the oo integration near 
lOq. Therefore, the next step toward evaluating the integral in (40) analytically is to expand 
the transmission coefficient T te (oou) in a power series about oo = oo 0 . Because our primary 
interest lies in evaluating (40) for a lossless slab with relative permittivity and permeability 
equal to —1 at the chosen frequency oo = oo 0 , that is, e(cu 0 )/e 0 = fi(u 0 )/fi 0 = —1, a power 
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series expansion of e(uj)/e 0 and which are needed in T te (uju), can be written as 


~~ — -1 H- — (a; 0 )(o; - u 0 ) + O [(uj - o> 0 ) 2 ] (41a) 

6q 6q CLUJ 

——- = -1 H- -^-{uq){u - u 0 ) + O [(a; - cn 0 ) 2 ] . (41b) 

fi 0 fi 0 auj 

For a passive material that is lossless at the frequency cc 0 (that is, e"(u 0 ) = fj”(uj 0 ) = 0), 
the frequency derivatives r jrj(u 0 ) and ^-(w 0 ) are also both zero because both e"(uj) and h"(lv) 
must be greater than or equal to zero in a passive material, and this is impossible in (41) for 
all uj near ca 0 unless ^y(^ 0 ) and ^-(w 0 ) are both zero. The coefficients J^o 0 )/e 0 an d 
of the (uj — u 0 ) terms in (41) are crucial to the evaluation of the power series for T te (uju) 
and the integral in (40). In fact, the final power series expansion for T te (uju ) obtained below 
shows that the time-domain solution differs substantially from the single-frequency solution 
only to the extent that ^(^ 0 )/ e o an d ^{^o)/Vo differ from zero. For lossless materials, both 
causality [20, sec. 84] and energy conservation [21, app. B] require that these coefficients 
have the following lower bounds 


1 de . .4 

(%) — j 
auj 


1 dfj, 4 

7 (^O) — 

dcu c<Jq 


(42) 


Consequently, the Jy(^ 0 )/ e o an d ^( U o)/Vo that vary the least rapidly near uj = uj 0 (and thus 
will reproduce the source fields most closely in the region z > 2 L) are given by 


k(uj) = = = _i + ±(uj-uj 0 ) + O[(uj- cc 0 ) 2 ] . (43) 

e o Vo ^0 

With (43) inserted into the expression for T te in (5), one finds that for the propagating 
waves (u 2 < /i 0 e 0 ) it can be approximated by 


T TE (cnn) « e 2lu ^ L , u 2 < ji 0 e 0 = 1/c 2 


(44) 


For the evanescent waves (u 2 > /i 0 e 0 ), the quantity k(uj)^ 0 (uu) / ^(uju) = —1 + 4[1 + 
l/(c 2 |Col 2 )](^ _ u o)/ UJ o + O [(a; ~ uj 0 ) 2 ] « -1 + 4(c v - uj 0 )/uj 0 + O [(a; - u 0 ) 2 ] if terms in 
1 / (c 2 1 Co 1 2 ) are neglected compared to unity and we find 


■^TE (k-’w) 


o^lCol-^ 


--4(o o — u 0 ) 2 + e -w|c 0 |n ’ 


u 2 > 1 /c 2 


(45) 
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Substitution of T te from (44) and (45) into (40) gives 


+oo 


E x (x,z,t) 


2vr 2 


du e 


—iujt 


sin [(a; — u; 0 )t 0 ] 


U — Un 


T(uy «[«®+c 0 (*- 2 i)] du 


u 2 <l/c 2 

sin [(a; — u 0 )t 0 \ T(u)e luJUX ~ u ^ z 


(uj-uj 0 ) e 2w|c ol L - A(u~u 0 ) 2 


du 


, z > d + L. (46) 


U 2 >1/c 2 

For large values of to, the delta-function approximation in (39) can be used to evaluate the 
frequency integration over the propagating spectrum in (46). Specifically, for large t 0 and 
\t±x/c\ < to (say \t±x/c\ < to/4), interchanging the order of the u and u integrations [27] 
reduces (46) to 


E x (x, z, t) 


p -iu 0 t r 

~ -- I r ^ e iu, 0 [ux+<: 0 (z-2L)} du 


2vr 


u 2 <l/c 2 

+oo 


2 vr 2 

u 2 >l /c 2 


T(u) 


sinhu; — u; 0 )t 0 ] e' lujux 

-ye f dudu, z>d + L. (47) 


(w-cu 0 ) e 2Hc o |i - ;4(w-^o) 2 


The remaining cu integration in (47) can be performed approximately by noting that for 
t 0 large and 1 1 — ux\ <C t 0 (say 1 1 — ux\ < to/4) [28], its major contribution comes from the 
integration of [sin(cu — u 0 )to]/(u — u 0 ) between u = u 0 — ir/to and u = u 0 + 7r/to, in which 
domain [sin(a; — u^)to\/{u — cu 0 ) can be replaced by its average value of approximately to/2. 
Then the u integration evaluates as 

+oo 

r sin[(o; — u 0 )t 0 ] e ' luJUX ~ UJ ^ z _ ii0t 


I = 


(u-uj 0 ) e 2uj tto\ L - A(u~ u 0 ) 2 


■re 


du 


uj 0 +ir/t 0 


t giwQttx-WolCob g-iwot 




uJ 0 -ir/t 0 [ 


e 2u, oKol L — -4(u; — cu 0 ) 2 


^O^O ^nUg-faJnlCnKg-^) } n 


1 — 1 V 0 g-^olCoT 

27T 


1 4- e^oKoT 
2n 


(48) 


The evanescent spectrum can be truncated at a value Co = -2* that makes the magnitude 
of the right-hand side of (48) about equal to its value at Co — 0. To find this value of Z t , 
first approximate the relevant part of the expression in (48) by 


g-^o Zt{z-L) 


_ ‘Vo p—LOnZtL 

27T 


1 + tZ^ e -u 0 ZtL 

2 7T 


(Vo -a , 0 Z t z 


71 


(49) 
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Then we want 

-i ^0^0 

1 — 2 ?r 

-| u Q to 

^271 

which can be solved to give 


^o e -“° ZtZ « In 

7r 


In 


1 - 


2 ?r 

CJq^O 


1 + 

CJq^O 


47T 

ca 0 to 


(50a) 


Z t « — In fa; t 0 > 27T. (50b) 

cUqZ \2nJ 

Since we are in the region z > d+ L, the value of z can be set equal to its minimum value, 
d + L, in (50b) to obtain 

z ‘ x ^hrL) Xn {^t) ■ (50c) 

Also, the function In 11 — ca 0 to e _ ^°^°l L /(27r) | has such a weak singularity with respect to 
variation in u as the value of its argument approaches zero that the u integration of (48) 
in (47) about this singularity is also negligible compared to the rest of the integration. 
Consequently, the u 2 > 1/c 2 integration in (47) over the evanescent waves can be truncated 
at | Co | = \Ju 2 — 1/c 2 ~ Z t and (47) reduces to 

E x {x, z, t ) «- / T{u)e iuJ ^ ux+ ^ z - 2L)] du , \t\ < t 0 /4, z > d + L. (51) 

27T J 
|c 0 l< 2 t 


The condition 1 1 — ux\ < t 0 /4 has been changed to \t\ < t 0 /A in (51) because for large t 0 the 
maximum value of ux is approximately equal to x \ Z t , which is much smaller than t 0 if t 0 
is large enough. 

If the time interval that the cos(u; 0 f) time dependence of the line source is turned on is 
shifted from [— t 0 ,t 0 \ to [0, 2t 0 ], the condition \t\ < to/4 in (51) shifts to 3t 0 /4 < t < 5t 0 /4. 
Moreover, since (51) would then hold for t > 3t 0 /4, causality demands that (51) would also 
hold if the signal turned off at 3fo/4 instead of 2f 0 . In other words, if the line source turned 
on at t = 0 with time dependence cos(o; 0 f), then t can replace 3f 0 /4 rs t 0 in (51) so that 

E x (x, z, t) ~ -- / T{u)e iuJ ° [ux+ ^ z ~ 2L)] d U , z > d + L (52) 

27r J 

ICol <Zt 


with 


Z t 


In I ^ 


ix 0 (d+L) \ 27r 


u n t S> 2n. 


(53) 


This equation can be recast in the form of (21) by returning the integration variable u to 
h/ujQ to get 


E x (x, z, t) 


g iio^t 

2 tt 


+H t 


T 0 {h)^ hx+ ^ z ~ 2L)] dh , z > d + L 


(54) 


-H t 
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with 


H t = 


U 0 Z t ) 2 + &ni 


00 


d + L 


1,11 ^ 

2ir 


+ k 2 m . 


(55) 


Here, 7 0 = (/cq 0 — h 2 )i and /cqo = u; o/ i o e o- I 11 th e region z > 2L, the value of z can be chosen 
in (50b) equal to its minimum value of 2 L in this region, so that H t in (55) is reduced in 
value to 

/Fi /mJ\1 2 

(56) 


H, 


■ In 


2 vr 


+ /cqq , z > 2L 


The field in (54) is merely that of the continuous-wave single-frequency held with its 
evanescent spectrum truncated at \h\ = H t given in (55) for d + L < z < 2L and (56) for 
z > 2L. Whereas the evanescent spectrum in (21) was truncated because of the presence 
of a small loss 5" in the slab, here in (54) the held of the lossless slab has its evanescent 
spectrum truncated because the sinusoidal time dependence of the line source turns on at 
time t — 0 instead of t — — oo and thus has had only a finite amount of time t to generate 
the evanescent waves. 

The resolution Arc just to the right of z = 2L is now a function of the time t that the 
source has been turned on. It is given from (56) as (see also (24)) 


Ax 


1.537T 

H t 


1.53t r 


[fi"(/o*)] 2 + £ 


(57) 


with a resolution enhancement of 


«4>o 


An 


27 tL 


In (fot) 


+ 1 


(58) 


where we have rewritten A; 0 o as 27 t/Ao and u 0 as 2nf 0 , Ao being the free-space wavelength 
and /o the cyclic frequency of the sinusoidal excitation. To attain a resolution enhancement 
R e just to the right of z = 2L, the line source would have to remain on for a time t given 
from (58) by 

1 „ L 7777 7 

(59) 


t 


l e 2 ^V^ 

fo 


This time is proportional to (and thus critically confirms) the estimate of the time obtained 
by Gomez-Santos [4] using a discrete split-frequency approximation (±Acn = u ± — cn 0 ) for a 
narrow band sinusoidal wave. 
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For example, to attain a resolution of R e = 5 at /o = 10 GHz, a line source illuminating 
an L — Ao wide slab would have to remain on for a time 

t ~ HT 10 e 2 ^' 724 « 39 minutes. (60) 

It is difficult to imagine a one wavelength wide slab made with material having small enough 
ohmic loss and structural inhomogeneities to maintain the build-up of evanescent fields for 
39 minutes at a frequency of 10 GHz. On the other hand, to attain a resolution of R e = 2.5, 
the same source would have to remain on for 

t « HT 10 e 27r 'V /(2 ' 5 ) 2-1 « 1.8 x HT 4 seconds. (61) 


These results indicate that for a one wavelength slab at a frequency of 10 GHz, a value of 
resolution enhancement R e around 2.5 may be feasible just to the right of z — 2 L. The 
resolution formulas (57)-(58) are confirmed by the numerical examples in Section IV. 

Comparing H t in (56) with H$ in (19), and (54) with (21), then allows us to determine 
an asymptotic approximation to (54) in the region d + L<z<2La,st —» oo by replacing 
1/8" in (31) with r = / 0 t to obtain for the y directed magnetic-current line source (and for 
z not too close to 2 L) 


E x (x,z,t ) 


r 2—z/L 


'x 1 + (2 L - z) 2 


In t — tan 


which confirms that the fields diverge to infinite values as t —» oo in the region d + L < z < 


2 L. 


B. Time-domain solution throughout the region to the right of the source (z > 0) 

The solution in (54) in the region z > d + L to the right of the lossless slab shows that 
the evanescent part of the transmitting spectrum of the slab is truncated by the value H t 
given in (55)-(56) that depends upon the amount of time (t) that the sinusoidal wave has 
been turned on (since the initial time t — 0). Therefore, for any finite time t, the fields of 
the lossless slab for all z > 0 to the right of the line source are given approximately by (8) 
with the infinite limits of integration (—oo, +oo) replaced by (—Ht, +H t ). This means that 
these fields of the lossless slab are finite everywhere for finite time t and yet approach the 
fields given in (10) as t —> oo. 
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If a very small loss is inserted into the slab material, an analysis similar to that performed 
in Section IIB shows that as t gets larger, the fields begin to approach the fields given in 
(10). As t gets even larger (eventually approaching infinity), however, the fields approach 
those given in (20) with the infinite fields in (20) replaced by large finite fields because 
of the very small but nonzero loss. After a long enough time t, the maximum resolution 
enhancement for the line source a distance d < L in front of the lossy slab is given just to 
the right of z = 2L by the formula in (25), that is 

*• “ \ Z ( 2^ W ") 2 + 1 (63) 

where 5" is the small loss defined as in (11) at the frequency u 0 . 

IV. NUMERICAL EXAMPLES 

The formulas for the resolution enhancement will be numerically validated in this section 
for both time-harmonic and time-domain line sources. Throughout, the central frequency 
is /o = 10 GHz and the slab width is L = Ao- We consider observation points only in the 
region z > 2 L, so the results are independent of the distance d between the line source 
and the slab as long as 0 < d < L. We begin with the expression (25) that determines the 
resolution enhancement for a lossy slab illuminated by a time-harmonic magnetic line source 
(zero bandwidth). The time-harmonic electric field from a single line source at the origin is 
obtained from (32) and (3) as 

+oo 

EJx, z) = / T TE (h) e ihx e iloz dh, z>d + L (64) 

27tA;oo J 

—oo 

where Tte {h) is given in (5) with the imaginary components of the permittivity and perme¬ 
ability satisfying (11). Figure 2 shows the magnitude of the integrand T^{h) e 2 * 7oi (corre¬ 
sponding to z — 2 L) as a function of h/koo for three different values of the loss parameter 5" 
(k 0 o is the free-space propagation constant evaluated at u> 0 ). The resolution enhancements 
obtained from (25) are shown as vertical lines and agree well with the observed effective 
spatial bandwidth of the integrand. 

Next consider a DNG slab that is lossless at the central frequency fo = 10 GHz and 
illuminated by a magnetic line source with time dependence v(t) = sin (cj 0 t) for 0 < t < T e 
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and v(t) = 0 otherwise. Hence the line source turns on at t — 0 and turns off at t = T e . We 
can then use (64) to get the following expression for the analytic-signal time-domain field 


E x (x,z,t) =— / W(h, z,t) e lhx dh, z>d+L 

2vr ./ 



( 66 ) 


(67) 


is the frequency spectrum corresponding to v(t). (The formula (65) can be obtained directly 
from (38) with in (36) replaced by V w in (67).) We can choose T e u 0 = 2nN where N is 
an integer to ensure that the ratio V u /u is bounded at uj = 0. However, in this numerical 
example the effective region of integration in ( 66 ) is confined to a narrow region centered 
on uo 0 that does not contain uj = 0. The time-domain spectrum ( 66 ) is a function of h that 
determines the spatial bandwidth of the time-domain field along a line perpendicular to the 


z axis. 

The calculation of the integral (66) is challenging because the integrand varies extremely 
rapidly near cu 0 . To ensure high accuracy we introduce a frequency-dependent loss that 
reduces the width of the region where the values of the integrand are non-negligiblc. This loss 
manifests itself in the (uj—lOq) 2 terms of the expressions for the permittivity and permeability 
(see (43)) 

= /fH = _! + + i 

do ^0 

which represents a perfectly lossless —1 DNG slab material at uj = lv 0 that satisfies the 
lower-bound requirements in (42). As noted earlier, a passive DNG material that is lossless 
at ui = lv 0 cannot have losses in the linear term of the power series expansion. 

The following nonuniform uj discretization is used to compute ( 66 ) for observation points 
along the line z = 2L + A 0 /1000. We let T e = 10 _3 s and evaluate the time-domain spectrum 
for 0 < h/koo < 3.5 using a 100000-point discretization. The step length varies as (uj — ca 0 ) 4 
for integration points near uj 0 . Away from uj 0 the step length is constant. The region of 
integration is 1 —10 -3 < uj/uj 0 < 1 + 10 -3 for 0 < h/k 00 < 2.5, and 1 —10 -9 < uj/uj 0 < 1 + 10 ~ 9 


1000(ca — u 0 ) \ 2 
/ 
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for 2.5 < h/koo < 3.5. This discretization may not be optimal, but it captures the variation 
of the integrand and ensures that the integral is computed accurately. 

Figure 3 shows the normalized magnitude of the analytic-signal time-domain spectrum 
| W(h, 2L + Ao/1000, t) | as a function of h/k 00 evaluated at three different times. Also plotted 
are the corresponding time-dependent resolution enhancements obtained from (58), which 
are seen to correctly predict the spatial bandwidth of the time-domain spectrum. 

We finally compute the analytic-signal time-domain electric field from two magnetic line 
sources that are A 0 /4 apart. Specifically, the line sources are located at (x,z) = (±A 0 /8,0). 
The constitutive parameters of the slab are given in (68) and the frequency spectrum for 
the time dependences of the line sources is given by (67). Figure 4 shows the normalized 
magnitude of the total analytic-signal electric field \E x (x,2L + A 0 /1000, t)\ along the line 
z = 2L + Ao/1000 at three different times. The positions of the line sources are indicated by 
two gray dots on the x axis of the figure. The resolution Ax predicted by (57) is Ax = 0.37A 0 , 
Ax = 0.32Ao, and Ax = 0.28Ao for t = 9 x 1 CT 6 s, t — 9 x 10 _5 s, and t — 9 x 10 ” 4 s, 
respectively. Since the actual distance between the line sources is Ax = 0.25Ao, we should 
expect to resolve the line sources only for the later time t — 9 x 1CT 4 s. Indeed, the electric 
field plot shows that the resolution improves with time and that the line sources are resolved 
only for t — 9 x 10 -4 s. For the two earlier times the line sources appear as a single peak. 
Thus, these numerical results confirm that the resolution enhancement increases with time 
and that the formula (57) gives a good estimate of the spatial resolution as a function of 
time. 
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FIG. 2: The magnitude of T^{h) e 2 * 7oL for a lossy slab as a function of h/k (M for 5" = 5.6 x 10 7 , 
5'" = 1.0 x 10~ 10 , and 5" = 4.3 x 10~ 14 . The corresponding values for the resolution enhancements 


are R e = 2.5, R e = 3.8, and R e = 5, as indicated with the vertical lines. 
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FIG. 3: The normalized magnitude of the analytic-signal time-domain spectrum as a function of 
h/koo f° r t = 10~ 6 s, t. = 10~ 5 s, and t = 10 _4 s. The corresponding values for the resolution 
enhancements are R e = 1.8, R e = 2.1, and R e = 2.5, as indicated with the vertical lines. 
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FIG. 4: The normalized magnitude of the analytic-signal time-domain electric field from two line 
sources as a function of x/Xo for t = 9 x 10 -6 s, t = 9 x 1CT 5 s, and t = 9 x 10“ 4 s, corresponding 
to a resolution of Ax = 0.37Ao, Ax = 0.32Ao, and Ax = 0.28Ao, respectively. The line sources are 
located at (x,z) = (±Aq/ 8,0) as indicated by the two gray dots. 
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V. CONCLUSION 


Plane-wave representations have been used to find the exact solution to a TE line source 
located in free space a distance d in front of a lossless or lossy magneto dielectric slab of 
width L for both single-frequency (zero bandwidth) sinusoidal excitation and for the time- 
domain (nonzero bandwidth) excitation of a sinusoidal wave that begins at time t = 0. In 
the special case of a lossless slab with relative permittivity and permeability equal to — 1 ( 
the —1 DNG slab) and d < L, the single-frequency source fields incident upon the slab are 
perfectly reproduced in the free-space region z > 2L to the right of the slab, but the fields 
to the right of the slab in the region d + L < z < 2L diverge to infinite values. These results 
are shown to hold regardless of whether the loss approaches zero before or after the limits 
of integration of the evanescent spectrum approach infinity. 

In contrast to the single-frequency fields, the time-domain (nonzero bandwidth) fields 
of the lossless —1 DNG slab for the sinusoidal excitation that turns on at t — 0 (rather 
than at t = —oo for the single-frequency excitation) remain finite everywhere for finite 
present time t and approach the fields of the single-frequency excitation only as t —> oo. 
In particular, perfect focusing (zero resolution of the line source) is not attainable after a 
finite time t because of the nonzero bandwidth [4], [ 6 ], [7] and the restrictions on the slowest 
possible variation with frequency of the permittivity and permeability allowed by causality 
and energy conservation. These results imply that the divergent infinite fields encountered 
in the lossless single-frequency solution to the —1 DNG slab are caused by the infinite energy 
in the single-frequency continuous-wave sinusoid that is imparted during the infinite amount 
of time between t = —oo and the present time t to the evanescent fields in the vicinity of 
the slab. Once the signal is made to turn on at some initial time t = 0, as would any signal 
in the laboratory, the fields remain finite everywhere for all future time t ^ oo, and thus, 
in principle, there appears to be no inconsistency inherent in postulating an ideal lossless 
infinitely long magneto dielectric slab with e(u; 0 )/e 0 = /j(u 0 )/fi 0 = —1. In reality, of course, 
finite losses, inhomogeneities in the material structure, noise levels, and the finite length of 
the slab, in addition to the nonzero bandwidth of the source coupled with the restrictions 
imposed by causality and energy conservation, will limit the resolution obtainable in the 
laboratory. 

The major frequency-domain and time-domain theoretical results were confirmed by the 
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direct numerical computations of Section IV. 
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by Ax ~ 1.28n/Hs ~ 1.287r/-^/(ln 5"/ L) 2 + (27t/Ao) 2 , which agrees well with the expression 
for the minimum H s y resolution obtained by replacing {e, /a}'' in [17, eq. (4)] with {J", J"} and 
solving for A (our Ax) to get Ax « 1.27r/Y / [ln((5 // /2)/L] 2 + (27r/Ao) 2 . 

[27] The order of integrations in (46) can be interchanged because (46) is absolutely integrable [22, 
p. 37], 

[28] Because the integration of the evanescent spectrum extends to u = Too, it appears that 
the inequality 1 1 — ux\ < to/4 cannot be satisfied for all u (unless x = 0). However, the 
following evaluation of the evanescent part of the integrals in (46) shows that truncating the 
u integration to finite limits introduces negligible error. 
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